Giant backscattering peak in angle-resolved Andreev reflection 
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It is shown analytically and by numerical simulation that the angular distribution of Andreev 
reflection by a disordered normal-metal-superconductor junction has a narrow peak at the angle 
of incidence. The peak is higher than the well-known coherent backscattering peak in the normal 
state, by a large factor G/Go (where G is the conductance of the junction and Go — 2e^/h). The 
enhanced backscattering can be detected by means of ballistic point contacts. 
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Coherent backscattering is a fundamental effect of 
time-reversal symmetry on the reflection of electrons by 
a disordered metal [0,1). The angular reflection dis- 
tribution has a narrow peak at the angle of incidence, 
due to the constructive interference of time-reversed se- 
quences of multiple scattering events. At zero tempera- 
ture, the peak is twice as high as the background. Coher- 
ent backscattering manifests itself in a transport experi- 
ment as a small negative correction of order Go = 2e^//i 
to the average conductance G of the metal {weak localiza- 
tion [^). Here we report the theoretical prediction, sup- 
ported by numerical simulations, of a giant enhancement 
of the backscattering peak if the normal metal (N) is in 
contact with a superconductor (S). At the NS interface 
an electron incident from N is reflected either as an elec- 
tron (normal reflection) or as a hole (Andreev reflection) . 
Both scattering processes contribute to the backscatter- 
ing peak. Normal reflection contributes a factor of two. 
In contrast, we find that Andreev reflection contributes 
a factor G/Gq, which is » 1. 

If the backscattering peak in an NS junction is so large, 
why has it not been noticed before in a transport experi- 
ment? The reason is a cancellation in the integrated an- 
gular reflection distribution which effectively eliminates 
the contribution from enhanced backscattering to the 
conductance of the NS junction. However, this cancel- 
lation does not occur if one uses a ballistic point contact 
to inject the current into the junction. We discuss two 
configurations, both of which show an excess conductance 
due to enhanced backscattering which is a factor G/Gq 
greater than the weak-localization correction. 

We consider a disordered normal-metal conductor 
(length L, width W, mean free path with N propa- 
gating transverse modes at the Fermi energy E-p) which 
is connected at one end to a superconductor (see inset of 
Fig. 1). An electron (energy Ep) incident from the oppo- 
site end in mode m is reflected into some other mode n, 
either as an electron or as a hole, with probability am- 
plitudes and rjjj^, respectively. The N x N matrices 
and r^'^ are given by 
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The Sij's are submatrices of the scattering matrix S of 
the disordered normal region. 
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where u, v, u' , v' are N x N unitary matrices, TZ = 1 — T, 
and T is a diagonal matrix with the transmission eigen- 
values Ti, T2, . . . Tm on the diagonal. 

We first consider zero magnetic field {B — 0). Time- 
reversal symmetry then requires that S* is a symmetric 
matrix, hence u' = , v' = . Eq. (1) simplifies to 
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We seek the average reflection probabilities ([r„„i|^), 
where (• • •) denotes an average over impurity configura- 
tions. Following Mello, Akkermans, and Shapiro [^, we 
assume that u is uniformly distributed over the unitary 
group. This "isotropy assumption" is an approximation 
which ignores the finite time scale of transverse diffusion. 
The refiection probabilities contain a product of four it's, 
which can be averaged by means of the formula [p[ 
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In the metallic regime N ^ L/l ^ I. In this large- 
N limit we may factorize (Efe=tfc' "^fc '''*;') i^^o {J2k'''k)^y 
which can be evaluated using M 

(EkfiTk)) - (A^//i)/rrf^ /(!/ cosh^ x). (5) 
The result for normal reflection is 

{K,n\')=ii+Snrn)N-\l^^l/L). (6) 

Off-diagonal (n ^ m) and diagonal {n = m) reflection 
differ by precisely a factor of two, just as in the normal 
state 1^. In contrast, for Andreev reflection we flnd 
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Off-diagonal and diagonal reflection now differ by an or- 
der of magnitude Nl/L ~ G/Gq > 1. 

Eqs. and (0) hold for B — Q. If time-reversal sym- 
metry is broken (by a magnetic field B > B^ = hJeLW) , 
then the matrices u, u' , v, v' are all independent ||7|. Car- 
rying out the average in the large- iV limit, we find 

(IC™P> = N-\l - ^l/L), (|r,tJ') = ^l/NL. (8) 

Diagonal and off-diagonal reflection now occur with the 
same probability. 

We have checked this theoretical prediction of a giant 
backscattering peak by a numerical simulation along the 
lines of Ref . ||9[] . The disordered normal region was mod- 
eled by a tight-binding Hamiltonian on a two-dimensional 
square lattice (dimensions 300 x 300, N = 126), with a 
random impurity potential at each site {L/l — 9.5). The 
scattering matrix S was computed numerically and then 
substituted into Eq. (1) to yield r°'^ and r^^. Results are 
shown in Fig. 1. This is raw data from a single sample. 
For normal reflection (bottom panel) the backscattering 
peak is not visible due to statistical fluctuations in the re- 
flection probabilities (speckle noise). The backscattering 
peak for Andreev reflection is much larger than the fluc- 
tuations and is clearly visible (top panel). A magnetic 
flux of lOh/e through the disordered region completely 
destroys the peak (middle panel). The arrow in the top 
panel indicates the ensemble-averaged peak height from 
Eq. (^, consistent with the simulation within the statis- 
tical fluctuations. The peak is just one mode wide, as 
predicted by Eq. (0) . UW>L the isotropy assumption 
breaks down Q and we expect the peak to broaden over 
W/L modes. Fig. 1 tells us that for L = W the isotropy 
assumption is still reasonably accurate in this problem. 

Coherent backscattering in the normal state is inti- 
mately related to the weak-localization correction to the 
average conductance. We have found that the backscat- 
tering peak for Andreev reflection is increased by a factor 
G/Gq. However, the weak- localization correction in an 
NS junction remains of order Go The reason is that 

the conductance 
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contains the sum over all Andreev reflection probabilities 
pof , so that the backscattering peak is averaged out. In- 
deed, Eqs. (0) and (p|) give the same G, up to corrections 
smaller by factors l/N and l/L. In order to observe the 
enhanced backscattering in a transport experiment one 
has to increase the sensitivity to Andreev reflection at the 
angle of incidence. This can be done by injecting the elec- 
trons through a ballistic point contact (width <C /, 
number of transmitted modes A'o). For B = 0, one can 
compute the average conductance from m 
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The density of transmission eigenvalues p{T) is known 
[|l|,0, in the regime Nq :S> 1, N ^ L/l. One finds 
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+ sind) = (NoL/Nl) cos^, i9e(0,7r/2). (lib) 

In the absence of time-reversal symmetry {B > Be) we 
find from the large- limit of Eqs. (1) and (H) that 
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This is just the classical addition in series of the Sharvin 
conductance NqGq of the point contact and the Drude 
conductance {NI/L)Gq of the disordered region. 

In Fig. 2 we have plotted the difference AG = {G{B = 
0)) - (G{B > Be)) of Eqs. (0) and If Nq/N < 

Z/L <C 1 the conductance drops from 2NoGo to NqGq 
upon breaking time-reversal symmetry. A doubling of the 
contact conductance at i? = is well-known in ballis- 
tic NS junctions {I ^ L). There it has a simple classical 
origin: An electron injected towards the NS interface is 
reflected back as a hole, doubling the current through 
the point contact. Here we find that the conductance 
doubling can survive multiple scattering by a disordered 
region (/ <C i), as a result of enhanced backscattering at 
the angle of incidence. 

As a second example we discuss how enhanced 
backscattering manifests itself when electrons are in- 
jected into a Josephson junction. The system considered 
is shown schematically in Fig. 3. A disordered metal 
grain is contacted by four ballistic point contacts (with 
Ni modes transmitted through contact i = 1,2, 3, 4). The 
scattering matrix S has submatrices Sij, the matrix el- 
ement Sij^nm being the scattering amplitude from mode 
m in contact j to mode n in contact i. The grain forms a 
Josephson junction in a superconducting ring. Coupling 
to the two superconducting banks is via point contacts 
3 and 4 (phase difference (j), same electrostatic poten- 
tial). Contacts 1 and 2 are connected to normal met- 
als (potential difference V). A current / is passed be- 
tween contacts 1 and 2 and one measures the conductance 
G = I/V as a function of (p. Spivak and Khniel'nitskii 
computed {G(6)) at temperatures higher than the Thou- 
less energy u^. They found a periodic modulation of 
the weak-localization correction, with amplitude of order 
Gq. Zaitsev and Kadigrobov et al. have discovered that 
at lower temperatures the amplitude increases to become 
much greater than Go ^6 17 1. Here we identify enhanced 
backscattering as the origin of this increase. 

The reflection matrices r'^° and r*^" (with elements 
fij,nm) contain the combined effect of scattering in the 
normal grain (described by the matrix S) and Andreev 
reflection at the two contacts with the superconductor. 
By summing a series of multiple Andreev reflections we 
obtain expressions analogous to Eq. (1), 



(10) 



(13a) 
(13b) 



2 



where we have abbreviated 



\ S21 S22 
( S31 



( Sl3 Si4 \ / S33 S34 \ 

\ S23 524^' V '^43 S44 ) ' 



S32 
S42 



e-'-^/^ 



The four-terminal generahzation of Eq. (0) is pi 
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Following Ref. we evaluate (G) by averaging S 

over the circular ensemble. At _B = this means that S — 
UU^ with U uniformly distributed in the group U{M) 
of M X M unitary matrices {M = Ei=i^i)- This is 
the circular orthogonal ensemble (COE). If time- reversal 
symmetry is broken, then S itself is uniformly distributed 
in U{M). This is the circular unitary ensemble (CUE). 
In the CUE we can do the average analytically for any 
Ni and 6. The result is 
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independent of 0. In the COE we can do the average 
analytically for 1 and = 0, and numerically |2y] 

for any Ni and 0. We find that the difference AG(0) — 
(G(0))coE - (G)cuE is positive for = 0, 
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with p = {N3 + N4)/{Ni + N2). The excess conductance 
( p^ ) is a factor G/Gq greater than the negative weak- 
localization correction, which is observable in Fig. 3 at 
(j) — TT. For Ni > 10 the finite-A^ curves (solid) are close 
to the large- limit (dotted) which we have obtained 
using the Green's function formulation of Refs. p3| , p^ . 

The excess conductance is a direct consequence of en- 
hanced backscattering. This is easiest to see for the 
symmetric case Ni = N2 = N, when (R^^) = 

{R\i) — (i?22)- Current conservation requires i?ii+-R2i + 
Rll + Rf^ ^N. For Af > 1 we may replace {f{Rtj)) by 
f{{Rij)). The average of Eq. ( |l^ then becomes 

(G/Go) = iA^ - - i?2i) + ^« - R2i)- (17) 

The first term ^A is the classical series conductance. 
The second term is the weak-localization correction due 
to enhanced backscattering for normal reflection. Since 
— i?|i) = 0(1) this negative correction to ^A^ can 
be neglected if A' ^ 1. The third term gives the excess 
conductance due to enhanced backscattering for Andreev 
reflection. Since {R^i — R21) = 0{N) this positive con- 
tribution is a factor G/Gq = 0{N) greater than the neg- 
ative weak-localization correction. 



In conclusion, we have predicted (and verified by nu- 
merical simulation) an order G/Gq enhancement of co- 
herent backscattering by a disordered metal connected 
to a superconductor. The enhancement can be observed 
as an excess conductance which is a factor G/Gq greater 
than the weak-localization correction, provided ballistic 
point contacts are used to inject the current into the junc- 
tion. The junction should be sufhciently small that phase 
coherence is maintained throughout. Several recent ex- 
periments p2| are close to this size regime, and might 
well be equipped with ballistic point contacts. 
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FIG. 1. Numerical simulation of a 300 x 300 tight-binding 
model for a disordered normal metal (L = 9.5/), in series 
with a superconductor (inset). The histograms give the modal 
distribution for reflection of an electron at normal incidence 
(mode number 1). The top two panels give the distribution 
of reflected holes (for B = and B = lO/i/eL^), the bottom 
panel of reflected electrons (for B — 0). The arrow indicates 
the ensemble-averaged height of the backscattering peak for 
Andreev reflection, predicted from Eq. (vn. 



FIG. 2. Excess conductance AG = {G{B = 0)) - 
{G{B > Be)) of a ballistic point contact in series with a dis- 
ordered NS junction (inset), computed from Eqs. (|l^) and 
(p^. At _B = the contact conductance is twice the Sharvin 
conductance NqGo, provided NqL/NI <^ 1. 

FIG. 3. Solid curves: excess conductance AG = (G((/'))coe 
— (G)cue of a four-terminal Josephson junction (inset), com- 
puted [|o) from Eqs. ^ and (|l^ for Ni = N2 = N, 
N3 = N4 = pN, with iV = 10. The dotted curves are the 
large- A'' limit The excess conductance &t 4> — Q is a, fac- 
tor G/Go = 0{N) larger than the negative weak-localization 
correction at cj) = n. 
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